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These lectures contain an elementary introduction to heavy quark symmetry and 
the heavy quark expansion. Applications such as the expansion of heavy meson 
decay constants and the treatment of inclusive and exclusive semilcptonic B decays 
are included. The use of heavy quark methods for the extraction of | V^j, | and \V u t\ 
is presented is some detail. 



1 Heavy Quark Symmetry 

In these lectures I will introduce the ideas of heavy quark symmetry and the 
heavy quark limit, which exploit the simplification of certain aspects of QCD 
for infinite quark mass, toq — > 00. We will see that while these ideas are 
extraordinarily simple from a physical point of view, they are of enormous 
practical utility in the study of the phenomenology of bottom and charmed 
hadrons. One reason for this is the existence not just of an interesting new 
limit of QCD, but of a systematic expansion about this limit. The technol- 
ogy of this expansion is the Heavy Quark Effective Theory (HQET), which 
allows one to use heavy quark symmetry to make accurate predictions of the 
properties and behavior of heavy hadrons in which the theoretical errors are 
under control. While the emphasis in these lectures will be on the physical 
picture of heavy hadrons which emerges in the heavy quark limit, it will be 
important to introduce enough of the formalism of the HQET to reveal the 
structure of the heavy quark expansion as a simultaneous expansion in powers 
of Aqcd/toq and a s (mg). However, what I hope to leave you with above all 
is an appreciation for the simplicity, elegance and coherence of the ideas which 
underlie the technical results which will be presented. The .interested reader 
is also encouraged to consult a number of excellent reviews p which typically 
cover in more detail the material in these lectures. 
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1.1 Introduction 



We begin by recalling the properties of charged current interactions in the 
Standard Model. They are mediated by the interactions with the bosons, 
which for the quarks take the form. 



(0 c t)r(l-j 5 )V C KM W^+h.c. (1) 

The 3x3 unitary matrix Vckm is 

Vckm = V cd V cs V cb I . (2) 



v ud 


V us 


v ub 


v cd 


Vcs 


Vcb 


Vtd 


Vts 


v tb 



The elements of Vckm have a hierarchical structure, getting smaller away from 
the diagonal: V u d, V cs and V t b are of order 1, V us and V c d are of order 10 _1 , 
V c b and V ts are of order 10 -2 , and V u b and Vtd are of order 10 -3 . By contrast, 
except for small effects associated with neutrino masses, the interaction of the 
W with the leptons is flavor diagonal. 

The CKM matrix is of fundamental importance, because it is the low en- 
ergy manifestation of the higher-energy physics which breaks the global flavor 
symmetries of the Standard Model. In the absence of the Yukawa couplings, 
the quark sector of the Standard Model may be characterized by its gauge sym- 
metry SU(3) x SU(2) x 17(1), and its global symmetry U(3)q x U(3)u x U(3) d 
which rotates the triplets of colored fields Q l L , U R and D R among each other. 
The global symmetry group has a total of 3 x 3 2 = 27 generators. With the 
addition of fields <f) and <p (possibly composite) carrying Higgs and conjugate- 
Higgs quantum numbers, one may write Yukawa-type interactions, 



A!? 



H (QituQ+Xj (Qi<PD j R ) +h.c. (3) 

which break the flavor symmetries explicitly. The complex matrices AJf and 
correspond to 2 x (2 x 3 2 ) = 36 independent parameters. They break the 
global flavor symmetries completely, except for a remaining conserved baryon 
number U(\)b- The 27 — 1 = 26 broken generators may be used to rotate 26 
of the Yukawa couplings to zero, leaving 36 — 26 = 10 physical parameters. 

When <j> and <f> get vacuum expectation values, one may go to the mass 
eigenstate basis for the quark fields, in which case the ten parameters are 
six quark masses and four parameters characterizing Vckm- An independent 
examination of Vckm confirms this counting. A 3 x 3 unitary matrix has 9 
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parameters, of which (2) = 3 are angles and the remaining 6 are complex 
phases. However, one may adjust 5 relative phases of the mass eigenstate 
quark fields, leaving 6 — 5 = 1 physical phase in Vckm- Thus Vckm is indeed 
characterized by four parameters, one of which is a CP-violating phase. 

It is an instructive exercise, left to the reader, to perform the analogous 
counting for the general case of U(Nj)q x U{Nj)u x U(Nf)jj global flavor 
symmetry. One finds that the Yukawa couplings contain Nf + 1 physical pa- 
rameters, of which 2Nf are quark masses. The remaining (Nf — l) 2 parameters 
characterize Vckm, with ^Nf(Nf — 1) being angles and \(Nf — l)(Nf — 2) be- 
ing complex phases. In particular, one notes that for Nf = 2 there is no CP 
violation in the weak interactions. 

Why is an understanding of QCD crucial to the study of the properties of 
Vckm? As an example, consider semileptonic b decay, b — > c£&, from which 
one would like to extract \V c b\. This process is mediated by a four-fermion 
operator, 

0bc = G^_^ {1 _^ )b _^ {1 _^ 5)L (4) 

The weak matrix element is easy to calculate at the quark level, 
GFV c b 

7f 



Ajuark = (c£v\ O bc \b) = j^-udjpe)^ (1 - J 5 )u b (p b ) Ui(pe)j^(l - ^)v v (p v ) 



(5) 

However, -4 qU ark is only relevant at very short distances; at longer distances, 
QCD confinement implies that free b and c quarks are not asymptotic states 
of the theory. Instead, nonperturbative QCD effects "dress" the quark level 
transition b — > civ to a hadronic transition, such as 

B -> DID or B -> D*£9 or ... (6) 

(In these lectures, we will use a convention in which a B meson contains a b 
quark, not a b antiquark.) The hadronic matrix element Ahadron depends on 
nonperturbative QCD as well as on GpV c b, and is difficult to calculate from 
first principles. To disentangle the weak interaction part of this complicated 
process requires us to develop some understanding of the strong interaction 
effects. 

There are a variety of methods by which one can do this. Perhaps the most 
popular, historically, has been use of various quark potential modelsB While 
these models are typically very predictive, they are based on uncontrolled 
assumptions and approximations, and it is virtually impossible to estimate the 
theoretical errors associated with their use. This is a serious defect if one builds 
such a model into the experimental extraction of a weak coupling constant such 
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as V c b, because the uncontrolled theoretical errors then infect the experimental 
result. 

These are issues which are important for the extraction of all the elements 
of Vckm- Let us pause now to review our current experimental knowledge of 
each of the magnitudes. The results are taken from the Particle Data BookH 
(The phases of the matrix elements must be extracted from CP violating asym- 
metries, as discussed elsewhere at this school.) 

We start with the submatrix describing mixing among the first two gener- 
ations. The parameter \V uc t\ is measured by studying the rates for neutron and 
nuclear j3 decay. Here the isospin symmetry of the strong interactions may be 
used to control the nonperturbative dynamics, since the operator c?7 M (l — 7 5 )it 
which mediates the decay is a partially conserved current associated with a 
generator of chiral SU(2)l x SU(2)r. The current data yield 

\V ud \ =0.9735 ±0.0008, (7) 

so \V u d\ is known at the level of 0.1%. The parameter \V US \ is measured simi- 
larly, via K — > tt£D£ and A — * ptvt. Here chiral SU (3)Z x SU (3)r must be used 
in the hadronic matrix elements, since a strange quark is involved; because the 
m s corrections are larger, \V US \ is only known to 1%: 

IK, | =0.2196 ±0.0023. (8) 

The Vckm elements involving the charm quark are not so well measured. 
One way to extract \V CS \ is to study the decay D — > Kl + vi. Unfortunately, 
there is no symmetry by which one can control the matrix element {K\s^{l — 
7 5 )c|D), since flavor SU(4) is badly broken. One is forced to resort to models 
for these matrix elements. The reported value is 

\V CS \ = 1.04 ±0.16, (9) 

but it must be said that this error estimate is not on very firm footing, and 
should probably be taken to be substantially larger. 

An alternative is to measure V cs from inclusive processes at higher energies. 
For example, one can study the branching fraction for W + — ► cs, which can 
be computed using perturbative QCD. The result of a preliminary analysis is 

\V CS \ = 1.00 ±0.13, (10) 

consistent with the model-dependent measurement. In this case, however, the 
error is largely experimental, and is unpolluted by hadronic physics. Similarly, 
one extracts \V c d\ from deep inelastic neutrino scattering, using the process 
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Vfi + d — > c + fi . This inclusive process may be computed perturbatively in 
QCD, leading to a result with accuracy at the level of 10%, 

\V cd \ = 0.224 ±0.016. (11) 

The elements of T^ckm involving the third generation are, for the most 
part, harder to measure accurately. The branching ratio for t — > bl + v can be 
analyzed perturbatively, but the experimental data are not very good. The 
present bound on \V t b\ is 

' F * fc ' -0.99 ±0.29. (12) 



\v td \ 2 + \v ts \ 2 + \v tb \> 

If one imposed the unitarity constraint \Vtd\ 2 + \Vt s \ 2 + \Vtb\ 2 — 1, then this 
would amount to a 15% measurement of |Vtb|, but this unitarity constraint is 
one of the properties of Vckm which one is trying to test. More generally, in 
fact, one should be wary of constraints on Vckm which impose unitarity as part 
of the analysis; while one often obtains tighter constraints in this way, these 
constraints have a different meaning than do direct measurements. What the 
comparison of a constrained and unconstrained "measurement" really tells you 
is whether the direct determination may be used to test the unitarity of Vckm • 

Unfortunately, there are as yet no direct extractions of \V tc i\ or \V ts \. One 
often speaks of these elements being measured in Bd — Bd and B s — B s mixing, 
but again, the hypothesis that the Standard Model is responsible for these 
processes is something that one really wants to check. The correct way to 
view this part of the experimental program is to say that the Standard Model, 
including the unitarity of Vckm, constrains Vt s and Vtd severely enough that 
testable predictions can be made for the mixing parameters Anid and Am s . 

This leaves us with the matrix elements V u b and V c b, for which we need an 
understanding of B meson decay. In these lectures we will discuss an approach 
to understanding the relevant hadronic physics which exploits the fact that the 
b and c quarks are heavy, by which we mean that mj, m c 3> Aqcd- The scale 
Aqcd is the typical energy at which QCD becomes nonperturbativc, and is of 
the order of hundreds of MeV. The physical quark masses are approximately 
rrib ~ 4.8 GeV and m c w 1.5 GcV. The formalism which we will develop will not 
make as many predictions as do potential models. However, the compensation 
will be that we will develop a systematic expansion in powers of Aqcd/™6,c, 
within which we will be able to do concrete error analysis. In particular, we 
will be able to estimate the error associated with the fact that m c may not be 
very close to the asymptotic limit m c >• Aqcd- Even where this error may be 
substantial, the fact that it is under control allows us to maintain predictive 
power in the theory. 
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1.2 The heavy quark limit 



Consider a hadron Hq composed of a heavy quark Q and "light degrees of 
freedom" , consisting of light quarks, light antiquarks and gluons, in the limit 
rriQ — > oo . The Compton wavelength of the heavy quark scales as the inverse of 
the heavy quark mass, Xq ~ 1/mg. The light degrees of freedom, by contrast, 
are characterized by momenta of order Aqcd, corresponding to wavelengths 
A^ ~ 1/Aqcd- Since A^ » Xq, the light degrees of freedom cannot resolve 
features of the heavy quark other than its conserved gauge quantum numbers. 
In particular, they cannot probe the actual value of Xq, that is, the value of 
m Q . 

We may draw the same conclusion in momentum space. The structure 
of the hadron Hq is determined by nonperturbative strong interactions. The 
asymptotic freedom of QCD implies that when quarks and gluons exchange 
momenta p much larger than Aqcd , the process is perturbative in the strong 
coupling constant a s (p) . On the other hand, the typical momenta exchanged 
by the light degrees of freedom with each other and with the heavy quark 
are of order Aqcd, for which a perturbative expansion is of no use. For these 
exchanges, however, p < ttlq, and the heavy quark Q does not recoil, remaining 
at rest in the rest frame of the hadron. In this limit, Q acts as a static source of 
electric and chromoelectric gauge field. The chromoelectric field, which holds 
Hq together, is nonperturbative in nature, but it is independent of ttlq. The 
result is that the properties of the light degrees of freedom depend only on the 
presence of the static gauge field, independent of the flavor and mass of the 
heavy quark carrying the gauge charge.^] 

There is an immediate implication for the spectroscopy of heavy hadrons. 
Since the interaction of the light degrees of freedom with the heavy quark is 
independent of toq, then so is the spectrum of their excitations. It is these 
excitations which determine the spectrum of heavy hadrons Hq. Hence the 
splittings Aj ~ Aqcd between the various hadrons Hq are independent of Q 
and, in the limit wiq — > oo, do not scale with uiq. 

For example, the bottom and charmed meson spectra are shown schemat- 
ically in Fig. 0, in the limit mj, m c ^> Aqcd- The light degrees of freedom are 
in exactly the same state in the mesons Bi and D^, for a given i. The offset 
Bi — Di = mi, — m c is just the difference between the heavy quark masses; in 
no way does the relationship between the spectra rely on an approximation 

We can enrich this picture by recalling that the heavy quarks and light 



"Top quarks decay too quickly for a static chromoelectric field to be established around 
them, so the simplifications discussed here are not relevant to them. 
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Figure 1: Schematic spectra of the bottom and charmed mesons in the limit mj,, m c S> Aqcd- 
The offset of the two spectra is not to scale; in reality, mj — m c ^> ~ Aqcd- 



degrees of freedom also carry spin. The heavy quark has spin quantum number 
Sq — which leads to a chromomagnetic moment 



HQ 



fJ 



2m 



Q 



(13) 



Note that /j,q — * as rriQ — * oo, and the interaction between the spin of the 
heavy quark and the light degrees of freedom is suppressed. Hence the light 
degrees of freedom are insensitive to Sq; their state is independent of whether 
Sq — ^ or Sq = Thus each of the energy levels in Fig. [I] is actually 

doubled, one state for each possible value of Sq. 

To summarize, what we see is that the light degrees of freedom are the 
same when combined with any of the following heavy quark states: 



l(T), QiU); Qa(T), QattY, 



(14) 



It is an 
A new 



where there are Nh heavy quarks (in the real world, Nh = 2) . T 
SU(2Nh) symmetry which applies to the light degrees of freedom 
symmetry means new nonperturbative relations between physical quantities 
It is these relations which we wish to understand and exploit. 
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The light degrees of freedom have total angular momentum Ji, which is 
integral for baryons and half-integral for mesons. When combined with the 
heavy quark spin Sq = |, we find physical hadron states with total angular 
momentum 

j=\j t ±h\. (is) 

If St 7^ 0, then these two states are degenerate. For example, the lightest heavy 
mesons have St = \ , leading to a doublet with J = and J = 1 . In the charm 
system we find that the states of lowest mass are the spin-0 D and the spin-1 
D*; the corresponding bottom mesons are the B and B*. The heavy quark 
spin operator Sq exchanges these two states. Writing the spin wave function 
\mQ,me), we have 

\M) = ^= (I Tl) - 1 IT)) , \M*(J 3 = 0)> = ±= (| Tl> + I IT)) • (16) 

Then it is easy to show that 

S Q \M) = ±\M*(J 3 = 0)), S Q \M*(J 3 = 0)) = \\M). (17) 

When effects of order 1 / toq are included, the chromomagnetic interactions 
split the states of given St but different J. This "hypcrfine" splitting is not 
calculable perturbatively, but it is proportional to the heavy quark magnetic 
moment /xq . This gives its scaling with vtiq : 

m D * — m D <~ l/m c 

rriB* — ins ^ 1/rrib . (18) 

From this fact we can construct a relation which is a nonperturbative prediction 
of heavy quark symmetry, 

m B » — m B = m 2 D , — m 2 D . (19) 

Experimentally, m B , — m B — 0.49 GeV 2 and m 2 D , — m 2 D = 0.55 GeV 2 , so 
this prediction works quite well. Note that this relation involves not just the 
heavy quark symmetry, but the systematic inclusion of the leading symmetry 
violating effects. 

Generally, the mass of a heavy hadron Hq may be expanded in inverse 
powers of tuq 

m(H Q )=m Q +K H + 0{l/m Q ), (20) 

where A B is independent of Q and is associated with the energy of the light 
degrees of freedom in the hadron Hq. For the lowest lying Jg — \ doublet, 
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this quantity is usually just referred to as A. From dimensional considerations, 
one expects A to be of order of a few hundred MeV. 

So far, we have formulated heavy quark symmetry for hadrons in their 
rest frame. Of course, we can easily boost to a frame in which the hadrons 
have arbitrary four-velocity = 7(1, v). For heavy quarks Q\ and Q2, the 
symmetry will then relate hadrons Hi(v) and H^iv) with the same velocity 
but with different momenta. This distinguishes heavy quark symmetry from 
ordinary symmetries of QCD, which relate states of the same momentum. To 
remind ourselves of this distinction, henceforth we will label heavy hadrons 
explicitly by their velocity: D(v), D*(v), B{v), B*{v), and so on. 

1.3 Semileptonic decay of a heavy quark 

Now let us return to the semileptonic weak decay b — > civ, but now consider 
it in the heavy quark limit for the b and c quarks. Suppose the decay occurs at 
time t = 0. For t < 0, the b quark is embedded in a hadron f^; for t > 0, the c 
quark is dressed by light degrees of freedom to H' c . Let us consider the lightest 
hadrons, il& = B(v) and H' c = D(v'). Note that since the leptons carry away 
energy and momentum, in general u/ti'. 

What happens to the light degrees of freedom when the heavy quark de- 
cays? For t < 0, they see the chromoelectric field of a point source with 
velocity v. At t = 0, this point source recoils instantaneously^] to velocity v'; 
the color neutral leptons do not interact with the light hadronic degrees of 
freedom as they fly off. The light quarks and gluons then must reassemble 
themselves about the recoiling color source. This nonperturbative process will 
generally involve the production of an excited state or of additional particles; 
the light degrees of freedom can exchange energy with the heavy quark, so 
there is no kinematic restriction on the excitations (of energy ~ Aqcd) which 
can be formed. There is also some chance that the light degrees of freedom 
will reassemble themselves back into a ground state D meson. The amplitude 
for this to happen is a function only of the inner product w = v • v' of the 
initial and final velocities of the color sources. This amplitude, £{w), is known 
as the Isgur-Wise functionta 

Clearly, the kinematic point v = v', or w = 1, is a special one. In this 
corner of phase space, where the leptons are emitted back to back, there is 
no recoil of the source of color field at t = 0. As far as the light degrees 
of freedom are concerned, nothing happens! Their state is unaffected by the 
decay of the heavy quark; they don't even notice it. Hence the amplitude for 
them to remain in the ground state is exactly unity. This is reflected in a 

''The weak decay occurs over a very short time St ~ l/Myy <S 1/Aqcd- 
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nonpcrturbative normalization of the Isgur-Wise function at zero recoil^ 

e(i) = i. (21) 

As we will see, this normalization condition is of enormous phenomenological 
use. It will be extremely important to understand the corrections to this result 
for finite heavy quark masses mb and, especially, m c . 

The weak decay b — > c is mediated by a left-handed current cj^(l — 7 5 )6. 
Not only does this operator carry momentum, but it can change the orientation 
of the spin Sq of the heavy quark during the decay. For a fixed light angular 
momentum Jg , the relative orientation of Sq determines whether the physical 
hadron in the final state is a D or a D* . However, the light degrees of freedom 
are insensitive to Sq , so the nonperturbative part of the transition is the same 
whether it is a D or a D* which is produced. Hence heavy quark symmetry 
implies relations between the hadronic matrix elements which describe the 
scmileptonic decays B — ► Dlv and B — > D*£p. 

It is conventional to parameterize these matrix elements by a set of scalar 
form factors. These arc defined separately for the vector and axial currents, 
as follows: 

(D(v')\c^b\B(v)) = h+{w)(v + v')> 1 + h_(w)(v -v'Y 
(D*(v',e)\c^b\B{v)) = h v (w)ie^ aP etv' a v fi 

{D(v')\&fT B b\B(v))=0 (22) 
(D*(v', e)| C7 M 7 5 6 \B(v)) = h Al (w)(w + l)e*" - e* ■ v[h Aa (w)v» + h A3 (w)v'»] . 

The set of form factors hi(w) is the one appropriate to the heavy quark limit. 
Other linear combinations are also found in the literature. In any case, the form 
factors are independent nonperturbative functions of the recoil or equivalently, 
for fixed mj, and m c , of the momentum transfer. However, in the heavy quark 
limit they correspond to a single transition of the light degrees of freedom, 
being distinguished from each other only by the relative orientation of the spin 
of the heavy quark. Hence they may all be written in terms of the single 
function £(w) which describes this nonperturbative transition. As we will 
derive later, the result is a set of relationsp 

h + (w) = h v (w) = h Al (w) = h Aa (w) = £(w) 

h-(w) =h A2 (w) =0, (23) 

which follow solely from the heavy quark symmetry. Of course, all of the form 
factors which do not vanish inherit the normalization condition (pi|) at zero 
recoil. This result is a powerful constraint on the structure of semileptonic 
decay in the heavy quark limit. 
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1.4 Heavy meson decay constant 

As a final example of the utility of the heavy quark limit, consider the coupling 
of the heavy meson field to the axial vector current. This is conventionally 
parameterized by of a decay constant; for example, for the B~ meson we 
define f B via 

(0\u^7 5 b\B-(p B ))=if B p B . (24) 

What is the dependence of the nonperturbative quantity f B on m B ? To ad- 
dress this question, we rewrite Eq. ( |24| ) in a form appropriate to taking the 
heavy quark limit, raj — > oo (which is equivalent to mi, —> oo). This entails 
making explicit the dependence of all quantities on m B . First, we trade the 
B~ momentum for its velocity, 

p B = m B v» . (25) 
Second, we replace the usual B~ state, whose normalization depends on m B , 
(B(p 1 )\B(p 2 ))=2E B 5^(p 1 ~p 2 ), (26) 
by a mass-independent state, 

\B(v)) = -±=\B(p B )), (27) 

satisfying 

{B{v 1 )\B{v 2 ))=2 1 5^{p 1 -p 2 ). (28) 
Then Eq. @ becomes 

y/mB~(0\ uj^ 5 b \B-(v)) = if B m B v^ . (29) 

The nonperturbative matrix element (0| u'y^j^b \B~(v)) is independent of m B 
in the heavy quark limit. Hence, we see that in this limit f B takes the form 

— 1/2 

f B = m B x (independent of ra B ) . (30) 

This makes explicit the scaling of f B with m B . It is more interesting to write 
this as axLcediction for the ratio of charmed and bottom meson decay constants. 
We findiffl 

^ = /^ +0 /'AqcD ) Aqcd\ 

jd V m s V m D m B J 

For the physical bottom and charm masses, of course, the correction terms 
proportional to Aqcd/wq could be important. 
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Figure 2: The nonleptonic decay of a b quark. 

2 Effective Field Theories 

2.1 General Considerations 

A central observation which underlies much of the theoretical study of B 
mesons is that physics at a wide variety of distance (or momentum) scales 
is typically relevant in a given process. At the same time, the physics at differ- 
ent scales must often be analyzed with different theoretical approaches. Hence 
it is crucial to have a tool which enables one to identify the physics at a given 
scale and to separate it out explicitly. Such a tool is the operator product ex- 
pansion, used in conjunction with the renormalization group. Here a general 
discussion of its application is given. 

Consider the Feynman diagram shown in Fig. ||[ in which a b quark de- 
cays nonleptonically. The virtual quarks and gauge bosons have virtualities fj, 
which vary widely, from Aqcd to My/ and higher. Roughly speaking, these 
virtualities can be classified into a variety of energy regimes: (i) /i 3> Mw] (ii) 
M w > fi » m b ; (iii) m b > fi > Aq CD ; (iv) fi « Aq CD . Each of these mo- 
menta corresponds to a different distance scale; by the uncertainty principle, 
a particle of virtuality \x can propagate a distance x ps 1/fi before being reab- 
sorbed. At a given resolution Ax, only some of these virtual particles can be 
distinguished, namely those that propagate a distance x > Ax. For example, 
if Ax > 1/Mw, then the virtual W cannot be seen, and the process whereby 
it is exchanged would appear as a point interaction. By the same token, as 
Ax increases toward 1/Aqcd, fewer and fewer of the virtual gluons can be 
seen explicitly. Finally, for \i sa Aqcd , it is inappropriate to speak of virtual 
gluons at all, because at such low momentum scales QCD becomes strongly 
interacting and an expansion in terms of individual gluons is inadequate. 

It is useful to organize the computation of a diagram such as is shown in 
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Fig. g in terms of the virtuality of the exchanged particles. This is important 
both conceptually and practically. First, it is often the case that a distinct 
set of approximations and approaches is useful at each distance scale, and one 
would like to be able to apply specific theoretical techniques at the scale at 
which they are appropriate. Second, Feynman diagrams in which two distinct 
scales yLii 3> [ii appear together can lead to logarithmic corrections of the 
form a™ ln™(/Lti/yU2), which for ln(/xi / ^2) ~ l/a; s can spoil the perturbative 
expansion. A proper separation of scales will include a resummation of such 
terms. 



2.2 Example I: Weak b Decays 

As an example, consider the weak decay of a 6 quark, b — ► cud, which is medi- 
ated by the decay of a virtual W boson. Viewed with resolution Ax < 1/Mw, 
the decay amplitude involves an explicit W propagator and is proportional to 



ef (I - 7 5 )6d 7/1 (l - j 5 ) u x <yl92 



p- 



(32) 



where p 11 is the momentum of the virtual W. Since nib *C Mw, the kinematics 
constrains p 2 <C M w , so the virtuality of the W is of order Mw, and it travels 
a distance of order 1/Mw before decaying. Viewed with a lower resolution, 
Ax > 1/Mw, the process b — > cud appears to be a local interaction, with 
four fermions interacting via a potential which is a 8 function where the four 
particles coincide. This can be seen by making a Taylor expansion of the 
amplitude in powers of p 2 /M w , 



C7 M (1 - 7 5 )6(i7 AI (l - 7 5 )w x 



M w 



1 



P 

M 2 
w 



p 



(33) 



The coefficient of the first term is just the usual Fermi decay constant, Gf/^/2. 
The higher order terms correspond to local operators of higher mass dimension. 
In the sense of a Taylor expansion, the momentum-dependent matrix element 
(|32|), which involves the propagation of a W boson between two spacetime 
points, is identical to the matrix clement of the following infinite sum of local 
operators: 



7 5 )6 



1 



(id) 2 + (id) 



d7/x(l-7 5 K 



(34) 



where the derivatives act on the entire current on the right. This expansion of 
the nonlocal product of currents in terms of local operators, sometimes known 
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as an operator product expansion, is valid so long as p 2 <C M w . For B decays, 
the external kinematics requires p 2 < m 2 , so this condition is well satisfied. 
In this regime, one may consider a nonrenormalizable effective field theory, 
with interactions of dimension six and above. The construction of such a low 
energy effective theory is also known as matching. As it is nonrenormalizable, 
the effective theory is defined (by construction) only up to a cutoff, in this case 
M\y- The cutoff is explicitly the mass of a particle which has been removed 
from the theory, or integrated out. If one considers processes in which one is 
restricted kinematically to momenta well below the cutoff, the nonrenormaliz- 
ability of the theory poses no technical problems. Although the coefficients of 
operators of dimension greater than six require counterterms in the effective 
theory (which may be unknown in strongly interacting theories), their matrix 
elements are suppressed by powers of p 2 /Af^. To a given order in p 2 /M w , 
the theory is well-defined and predictive. 

From a modern point of view, in fact, such nonrenormalizable effective 
theories are actually preferable to renormalizable theories, because the non- 
renormalizable terms contain information about the energy scale at which the 
theory ceases to apply. By contrast, renormalizable theories contain no such 
explicit clues about their region of validity. 

In principle, it is possible to include effects beyond leading order inp 2 /M w 
in the effective theory, but in practice, this is usually quite complicated and 
rarely worth the effort. Almost always, the operator product expansion is trun- 
cated at dimension six, leaving only the four-fermion contact term. Corrections 
to this approximation are of order m 2 /M w ~ 10~ 3 . 

2.3 Radiative Corrections 

At tree level, the effective theory is constructed simply by integrating out the 
W boson, because this is the only particle in a tree level diagram which is 
off-shell by order M w . When radiative corrections are included, gluons and 
light quarks can also be off-shell by this order. Consider the one-loop diagram 
shown in Fig. |^. The components of the loop momentum fc M are allowed to 
take all values in the loop integral. However, the integrand is cut off both in 
the ultraviolet and in the infrared. For k > M\y, it scales as d 4 k/k 6 , which 
is convergent as k — > oo. For k < nib, it scales as d 4 fc/fc 3 mf,M^, which is 
convergent as k — > 0. In between, all momenta in the range mt < k < M\y 
contribute to the integral with roughly equivalent weight. 

As a consequence, there is potentially a radiative correction proportional 
to a s ln(Mn// m fc)- Even if a s (/i) is evaluated at the high scale fi — Mw, such 
a term is not small in the limit Mw — > oo. At n loops, there is potentially a 
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Figure 3: The nonleptonic decay of a b quark at one loop. 



term of order a™ ln n (Afty/m&). For a s hx^Mw/mb) ~ 1, these terms need to 
be resummed for the perturbation series to be predictive. The technique for 
performing such a resummation is the renormalization group. 

The renormalization group exploits the fact that in the effective theory, 
operators such as 



receive radiative corrections and must be subtracted and renormalized. (Here 
the color indices i and j are explicit.) In dimensional regularization, this means 
that they acquire, in general, a dependence on the renormalization scale \i. Be- 
cause physical predictions are independent of /x, in the renormalized effective 
theory the operators must be multiplied by coefficients with a compensating 
dependence on /x. It is also possible for operators to mix under renormal- 
ization, so the set of operators induced at tree level may be enlarged once 
radiative corrections are included. In the present example, a second operator 
with different color structure, 



is induced at one loop. The interaction Hamiltonian of the effective theory is 
then 



(35) 



Oh = c j7 ^(1 - 7 5 )6^ d,7 M (l - 7 5 K , 



(36) 



n eS = C7(m)Oj( M ) + C n (ji)On(n) , 



(37) 



and it satisfies the differential equation 



(38) 
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By computing the dependence on \x of the operators Oi(/i), one can deduce 
the ^-dependence of the Wilson coefficients Ci(p). In this case, a simple 
calculation yields 



ChJ ^6/23 / , A/r nx -12/23' 

a s (M w )\ f a s (M w ^ 



(fi) J V "s(M) 



(39) 



For /j, — nib, these expressions resum all large logarithms proportional to 
aj In™ (M w /m b ). 

The decays which are observed involve physical hadrons, not asymptotic 
quark states. For example, this nonleptonic b decay can be realized in the 
channels B — > Dir, B — > D*tttt, and so on. The computation of partial decay 
rates for such processes requires the analysis of hadronic matrix elements such 
as 

(Dir\ c 7 '<(l - 7 5 )6w 7ai (1 - 7 5 )d \b) . (40) 

Such matrix elements involve nonperturbative QCD and are extremely difficult 
to compute from first principles. However, they have no intrinsic dependence 
on large mass scales such as Mw- Because of this, they should naturally be 
evaluated at a renormalization scale /i <C Mw, in which case large logarithms 
hi(M\nr/mb) will not arise in the matrix elements. By choosing such a low 



scale in the effective theory (37), all such terms are resummed into the coeffi- 
cient functions Ci(mb). As promised, the physics at scales near Mw has been 
separated from the physics at scales near nib, with the renormalization group 
used to resum the large logarithms which connect them. In fact, as we will 
see in the next section, nonperturbative hadronic matrix elements are usually 
evaluated at an even lower scale \i m Aqcd -C nib, explicitly resumming all 
perturbative QCD corrections. 



3 Heavy Quark Effective Theory 

We have already extracted quite a bit of nontrivial information from the heavy 
quark limit. We have found the scaling of various quantities with niQ, we have 
studied the implications for heavy hadron spectroscopy, and we have found 
nonperturbative relations among the hadronic form factors which describe 
semileptonic b — > c decay. However, all of these results have been obtained 
in the strict limit niQ — > oo. If the heavy quark limit is to be of more than 
academic interest, and is to provide the basis for quantitative phenomenology, 
we have to understand how to include corrections systematically. There are 
actually two types of corrections which we would like to include. Power correc- 
tions are subleading terms in the expansion in Aqcd / mQ\ those proportional to 
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Aqcd/w c are the most worrisome, because of the relatively small charm quark 
mass. Logarithmic corrections arise from the implicit dependence of quantities 
on rriQ through the strong coupling constant a s (rriQ) ~ 1/ I^toq/Aqcd)- For 
the physical values of mb and m c , either of these could be important. What 
we need is a formalism which can accommodate them both. 

In short, we need to go from a set of heavy quark symmetry predictions 
in the wiq — > oo limit, to a reformulation of QCD which provides a controlled 
expansion about this limit. The formalism which does the job is the Heavy 
Quark Effective Theory, or the HQET. The purpose of the HQET is to allow us 
to extract, explicitly and systematically, all dependence of physical quantities 
on toq, in the limit tuq 3> Aqcd- In these lectures, we will develop only enough 
of the technology to treat the dominant leading effects, providing indications 
along the way of how one would carry the expansion further. 

The HQET, as foim^atoEherje, was-dgveloped in a series of papers going 
back to the late 1980's^>Byy>BE^My which the reader who is interested 
in tracking its historical development may consult. 

3.1 The effective Lagrangian 

Consider the kinematics of a heavy quark Q, bound in a hadron with light 
degrees of freedom to make a color singlet state. The small momenta which Q 
typically exchanges with the rest of the hadron are of order Aqcd ttlq , and 
they never take Q far from its mass shell, p 2 q = voq. Hence the momentum pg 
can be decomposed into two parts, 

p£ = m Q v^ + k» , (41) 

where rriQV^ is the constant on-shell part of Pq, and fc M ~ Aqcd is the small, 
fluctuating "residual momentum" . The on-shell condition for the heavy quark 
then becomes 

m 2 Q = (mqv^ + k^) 2 = m Q + 2m Q v ■ k + k 2 . (42) 

In the heavy quark limit we may neglect the last term compared to the second, 
and we have the simple condition 

v ■ k = (43) 

for an on-shell heavy quark. Here the velocity v 1 * functions as a label; since 
soft interactions cannot change i> M , there is a velocity supers election rule in the 
heavy quark limit, and is a good quantum number of the QCD Hamiltonian. 
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We find the same result by taking the toq — > oo limit of the heavy quark 
propagator, 

i y l±± * (44) 

P — rriQ + it 2 v ■ k + ie 

In this limit the propagator is independent of tuq. The projection operators 

P ± = if* (45) 

project onto the positive (P+) and negative (P_) frequency parts of the Dirac 
field Q. This is clear in the Dirac representation in the rest frame, in which P + 
and P- project, respectively, onto the upper two and lower two components of 
the heavy quark spinor. In the limit tuq — > oo, in which Q remains almost on 
shell, only the "large" upper components of the field Q propagate; mixing via 
zitterbewegung with the "small" lower components is suppressed by l/2mg. 
Hence the action of the projectors on Q is 

P+Q(x) = Q(x) + 0{l/m Q ), P-Q(x)=0 + O(l/m Q ). (46) 

More precisely, these relations should be understood as pertaining to those 
modes of the field Q{x) which annihilate heavy quarks and antiquarks in a 
heavy meson. 

The momentum dependence of the field Q is given by its action on a heavy 
quark state, 

Q(x) \Q(p)) = e-** |0) . (47) 

If we now multiply both sides by a phase corresponding to the on-shell mo- 
mentum, 

j m * v - x Q{x) \Q( P )) = e~ ik - x |0) , (48) 

the right side of this equation is independent of toq. Hence the left side must 
be, as well. Combining this observation with the argument of the previous 
paragraph, we are motivated to define a mQ-independent effective heavy quark 
field h v (x), 

h v (x) =e imQV - x P+Q(x). (49) 

Note that the effective field carries a velocity label v and is a two-component 
object. The modifications to the ordinary field Q{x) project out the positive 
frequency part and ensure that states annihilated by h v (x) have no dependence 
on tuq. Hence, these are reasonable candidate fields to carry representations 
of the heavy quark symmetry. Of course, the small components cannot be 
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neglected when effects of order l/m,Q arc included. In the HQET they are 
represented by a field 

H v {x) =e imQV - x P-Q{x). (50) 

The field H v {x) vanishes in the vhq — ► oo limit. 

The ordinary QCD Lagrange density for a field Q{x) is given by 

£qcd = Q{x){ip-m Q )Q(x), (51) 

where = — igA^T a is the gauge covariant derivative. To find the La- 
grangian of the HQET, we substitute 

Q(x)^e- tm ^- x h v (x) + ... (52) 

into £qcd and expand. With the aid of the projection identity P + 7 A 'P + = , 
we find 

£hqet = h v (x)iv ■ Dh(x) . (53) 
This simple Lagrangian leads to the propagator we derived earlier, 

(54) 



v ■ k + it ' 

and to an equally simple quark-gluon vertex, 

igTVA; . (55) 

Note that both the propagator and the vertex are independent of mg, re- 
flecting the heavy quark flavor symmetry. They also have no Dirac structure, 
reflecting the heavy quark spin symmetry. Our intuitive statements about the 
structure of heavy hadrons have been promoted to explicit symmetries of the 
QCD Lagrangian in the limit rriQ — ► oo. 

It is straightforward to include power corrections to £hqet- Write Q(x) 
in terms of the effective fields, 

Q(x) = e~ im ^- x [h v (x) + H v (x)} , (56) 
and apply the classical equation of motion (ip — niQ)Q(x) = 0: 

iph v (x) + (ip - 2m Q )H v (x) = 0. (57) 
Multiplying by P_ and commuting f to the right, we find 

{iv ■ D + 2m Q )H v (x) =ip±h v (x), (58) 
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where Dj_ = — v^v ■ D. We then substitute Q(x) into £qcd as before, 
eliminate H v (x) and expand in 1/toq to obtain 

Ajqet = h v iv ■ Dh v + h v ip± 1 ip±h v 

iv ■ D + 2mg 



h v iv ■ Dh v + 



[ h v (iD ± ) 2 h v + | h v a aP G afi h v \ + ... . (59) 
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The leading corrections have a simple interpretation, which becomes clear in 
the rest frame, w M = (1,0,0,0). The spin-independent term is 

7^— Ok = ^— h v (iD ± ) 2 h v -» -— !— h v {iD) 2 h v , (60) 
2m Q 2m Q 2m Q 

which is the negative of the nonrelativistic kinetic energy of the heavy quark. 
Because of the explicit factor of 1/2toq, this term violates the heavy flavor 
symmetry. The spin-dependent part is 

^ Oa = JL f h v ^G atihv - JL h v a ij T a h v x gG% = g^B", (61) 

which is the coupling of the spin of the heavy quark to the chromomagnetic 
field. Because it has a nontrivial Dirac structure, this term violates both the 
heavy flavor symmetry and the heavy spin symmetry. For example, Oq is 
responsible for the D — D* and B — B* mass splittings. These correction terms 
will be treated as part of the interaction Lagrangian, even though Ok has a 
piece which is a pure bilinear in the heavy quark field. 

3.2 Effective currents and states 

The expansion of the weak interaction current C7 M (1 — j 5 )b is analogous. How- 
ever, here we must introduce separate effective fields for the charm and bottom 
quarks, each with its own velocity: 

b^h b v , c->h c v ,. (62) 

Then a general flavor-changing current becomes, to leading order, 

cTb^h c v ,Th b v , (63) 

where T is a fixed Dirac structure. With the leading power corrections, this is 
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Q i Q ? q q Q q 



Figure 4: Tree level plus the one loop renormalization of the current qTQ in QCD. The box 
represents the current insertion. 

The effective currents, and other operators which appear in the HQET, may 
often be simplified by use of the classical equation of motion, 

iv-Dh v (x)=0. (65) 

However, it is only safe to apply these equations naively at order I/tuq; at 
higher order the application of the equations of motion involves additional 
subtletiesBffl 

To complete the effective theory, we need mQ-independent hadron states 
which are created and annihilated by currents containing the effective fields. 
For example, there is an effective pseudoscalar meson state \M(v)) which cou- 
ples to the effective axial current qj^^hy, with a coupling Fm which is inde- 
pendent of m,Q , 

(0| q^~f 5 h v \M(v)) = iF M v» . (66) 
At lowest order, Fm is related to a conventional decay constant such as fg by 

fB = F M /VmB, (67) 
from which we immediately find the relationship ( |3l| ) between fu and fs- 

3.3 Radiative corrections 

We can use the effective Lagrangian £hqet to compute the radiative correc- 
tions to the matrix element d66|). In particular, we would like to extract the 
dependence of Fm on lnmg. This dependence comes through the one-loop 
renormalization of the current qj^j 5 Q. At lowest order, of course, the renor- 
malization is straightforward: we simply compute the set of graphs found in 
Fig. ^. The result is finite, because the current is (partially) conserved, and 
we extract a result of the form 

g 7 M 7 5 (l + 7o Hm Q /m q ) + ...). (68) 
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h v q h v q h r q h v q 

Figure 5: Tree level plus the one loop renormalization of the current q7 M 7 5 h l , in HQET. 
The double line is the propagator of the effective field h v . 

Note that there is no explicit dependence on the renormalization scale /i, since 
there is no divergence to be subtracted. 

The same result may be obtained in the effective theory. In this case we 
must match the currents in full QCD onto HQET currents of the form qj^j 5 ^. 
This step will induce a matching coefficient containing the explicit dependence 
on rriQ, which is absent, by construction, from the operators and Lagrangian of 
the HQET J] In addition, the effective current will not necessarily be conserved, 
since the ultraviolet properties of QCD and the HQET differ. Hence the form 
of the matching, once radiative corrections are included, is 

qYl 5 Q -> C(m Q , M ) x qr~f 5 hv(») ■ (69) 

We can deduce the form of C {uiq , /i) by considering the renormalization 
of the effective current 97 M 7 5 /i„, shown in the last three terms of Fig. g. These 
diagrams, computed in the effective theory, are independent of tuq. However, 
in general they are divergent, so they depend on the renormalization scale /i; 
the renormalization takes the form 

q^ 5 h v ((i) = q^ 5 h v (m q ) x (l +7o^ Mp/m q ) + ...), (70) 

where here m q acts as an infrared cutoff. The [i dependence in the second term 
must be canceled by C(rriQ, /i). Since the logarithm depends on a dimcnsionless 
ratio, C(mQ,/i) must be of the form 

C(m Q ,n) = 1 + 7o^; hi(m Q //i) + . . . . (71) 

Comparing the dependence on toq of C(rriQ,fi) and the expansion (p8|), we see 
that -yo = 70 ■ 

c In general, the matching procedure at order a s can also induce new Dirac structures qTh v . 
They do not affect the leading logarithms discussed here. 
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However, the effective theory allows us to go beyond leading order, and to 
resum all corrections of the form a™ In™ m,Q. We do this with the renormaliza- 
tion group equations, which express the independence of physical observables 
on the renormalization scale y,. In this case, they require that the [i depen- 
dence of C(m,Q, n) cancel that of the one loop diagrams in Fig. ^, under small 
changes in /i: 

^C7(m QlM ) = -70^. (72) 

The logarithms are resummed because the partial derivative is promoted to a 
total derivative with respect to /i, including the implicit dependence on /i of 
the coupling constant a s (/x): 

d d n . . d 

m = -^ + - 

2 25 

fa = 11 - -N f = - foriV / =4, (73) 

where Nf is the number of light flavors. We compute the anomalous dimension 
7o from the ultraviolet divergent parts of the one loop diagrams shown in Fig. ^. 

It is instructive to perform the radiative correction to the current in detail, 
since this is different from the diagrams one is used to in ordinary QCD. With 
the HQET Feynman rules, the diagram may be written in Feynman gauge as 

Cf fa) V / -^X v aa 7 V — v a u h x 4 > (74) 
J (2tt) 4 e i vq q A 

where Cf — | is the color factor. This expression may be simplified to 

f d 4 ~ e o q 13 

- C s ig 2 vt^ ^j 5 u h ^ / - , (75) 



(27r) 4 - e q 4 v ■ q 



which by Lorentz invariance is simply 



-C,iifrtf>fu h v< J (76) 

Rotating to Euclidean space, performing the integral and extracting the pole 
in e, we find 

va^ 5 u h x (2Cf) + finite . (77) 

lD7r z e 
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Since /i c = 1 + eln/i + . . ., the one loop contribution to the matrix element 
then depends on In jj, as 

v a ^ h u h x(2C f )^ In fi. (78) 

The contribution to the anomalous dimension is then 2Cf = f . 

The calculation of the wavefunction renormalization of the heavy quark 
(the third diagram in Fig. |^) is similar, but requires a novel version of the 
Feynman trick, 

1 '°° ^ (79) 



ab J (a + \b) 2 

One also has to pick out the term which cancels the l/v ■ k pole in the heavy 
quark propagator. With the factor of \ which accompanies the contributions 
from wavefunction renormaliaztion, the result is 

\ x vrf<f u h x(4C/)^ In /x. (80) 

A 47T 

Including the usual QCD renormalization of the light quark field, we find from 
the three terms in Fig. ^, respectivelypfl 

8 1 / 8\ 1 /16\ „ 
7o = x+o ho +o hr = 4- (81) 



3 2 

The solution of the renormalization group equation is 

c{mQ ^-\-^w) -{-^w) • (82) 

This then yields the leading logarithmic correction to the ratio Jb/}d'- 

6/25 



Ib _ /tod / a s {m c ) 



Jd V m B \a e (mb) 



(83) 



The radiative correction is approximately a ten percent effect. In fact, it has 
a simple physical interpretation. For virtual gluons of "intermediate" energy, 
m c < Eg < mj), the bottom quark is heavy but the charm quark is light. Such 
gluons contribute to the difference between /g and fo even in the heavy quark 
limit. 

In summary, then, the purpose of the HQET is to make explicit all depen- 
dence of observable quantities on tuq. The logarithmic dependence, through 
cKs(toq), arises from intermediate virtual gluons with m c < E g < nib- We 
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obtain these corrections by computing perturbatively with the HQET La- 
grangian, then using the renormalization group to resum the logarithms to 
all orders. The power dependence, I/ttoq, is extracted systematically in the 
heavy quark expansion. We have seen how to expand the Lagrangian and the 
states to subleading order; the application of the expansion to a physical decay 
rate will be presented in the next section. 

These lectures are meant to be pedagogical, so we will only treat the lead- 
ing corrections to a few processes. However, the state of the art goes signifi- 
cantly beyond what will be presented here. For many quantities, not only the 
leading logarithms, a™ In™ tuq, but the subleading (two loop) logarithms, of or- 
der a ra+1 In™ tuq, have been resummed. Similarly, many power corrections are 
known to relative order l/m,Q. It is particularly important phenomenologically 
to take into account the corrections of order l/m^. 

4 Exclusive B Decays 

We now have the tools we need for an HQET treatment of the exclusive semilep- 
tonic transitions B — > D lv and B — > D* lv. Earlier, we argued on physical 
grounds that in the heavy quark limit all of the hadronic matrix elements 
which appear in these decays are related to a single nonpcrturbative function 
£,{w). Now we will sharpen this analysis to actually derive these relations, and 
to include radiative and power corrections. In fact, almost all of our effort will 
go into the power corrections, since the radiative corrections to the transition 
currents are computed just as in the previous section. 

4-1 Matrix element relations at leading order 

The transitions in question require the nonperturbative matrix elements 

(D(v')\ c^b \B(v)) , (D*(v', e)\ c^b \B(v)) , (D*(v', e)| c 7 " 7 5 fc \B(v)) , 

(84) 

parameterized in terms of form factors as in Eq. (|22|) . Our first task is to 
derive the relations between these form factors, as promised earlier. These 
relations depend on the heavy quark symmetry, that is, on the fact that the 
spin quantum numbers of Q and of the light degrees of freedom are separately 
conserved by the soft physics. Hence we need a representation of the heavy 
meson states in which they have well defined transformations separately under 
the angular momentum operators Sq and J(. In particular, the representation 
must reflect the fact that a rotation by Sq can exchange the pseudoscalar 
meson M(v) with the vector meson M*(v, e). 
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Thpaolution is to introduce a "superfield" M{v), defined as the 4x4 Dirac 
matrix BE3 

M(v) = ±±t [rM;(v,e)--y 5 M(v)] = V(v,e) + P(v). (85) 

Under heavy quark spin rotations Sq, M.(y) transforms as 

M(v) -> D(S Q )M(v) , (86) 

and under Lorentz rotations A, as 

M(v) D{A)M(A- 1 v)D- 1 (A). (87) 

Here D(- ■ •) is the spinor representation of 50(3, 1). The superfield satisfies 
the matrix identity 

P+M(v)P- = M(v), (88) 

so it transforms the same way as the product of spinors h v q, representing a 
heavy quark and a light antiquark moving together at velocity v^. 

It is straightforward to verify the transformation properties of the super- 
field in the rest frame, in which i> M = (1, 0, 0, 0). In this frame, the spinor repre- 
sentation of the angular momentum operator J has components S l — 57 5 7°7 J - 
It acts on the superfield by 3 l M = [S l ,Ai]. Defining the polarization vectors 
<4 = (0, i/\/2,±i/V2,0) and e£ = (0,0,0, 1), it is easy to check that 

J 2 P = J 3 P = 0, J 2 V(e) = 2V(e) , J 3 V{e±) = ±V(e±) , J 3 y(e) = 0, 

(89) 

so P has spin zero and V(e) spin one. On the other hand the heavy quark spin 
operator Sq has the same component representation but acts only on the left, 
(SqYM — S*M. One may then check that 

(S Q ) 3 P = iu( £3 ), (S Q ) 3 U( £3 ) = ip. (90) 

As promised, heavy quark spin tranformations exchange the pseudoscalar and 
vector mesons. 

A current which mediates the decay of one heavy quark (Q) into another 
(Q') is of the form h V 'Th v . Under a rotation by Sq, the effective field h v 
transforms as 

h v ^D(S Q )h v , (91) 

while h v i is unchanged. The current would remain invariant if we took T to 
transform as 

T^TD-^Sq). (92) 
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On the other hand, the matrix element of superfields 



(M'(v')\h v ,Th v \M(v)) (93) 

is invariant if we rotate both h v and M.(y) by the same Sq. With the trans- 
formation law ( |86| ) for M(v), it follows that the S^-invariant matrix element 
must be proportional to TM. (v) . When we also consider rotations under Sq> , 
we find that the matrix element is restricted to the general form 

(M'(v')\ h v . T h v \M(v)) = -s/M M M w Tr \m'(v') FM(v) F («,t/)l ■ (94) 

The product of masses in front is a convention which restores the relativistic 
normalization of the states. Note that the heavy quark symmetry allows an 
arbitrary 4x4 Dirac matrix F(v,v') to act on the "light quark" part of the 
superfields. Its presence reflects the fact that only Lorentz symmetry constrains 
the spin of the light degrees of freedom during the decay. 

A general expansion of F(v, v') in terms of scalar functions Fi(w — v ■ v') 
takes the form 

F(v,v') =F 1 (w) + F 2 {w)f + F 3 (w)i>' + F A (w)W . (95) 

However, the identity (|8|) applied to the matrix element (|94|) yields 

F(v, v') = P- F(v,v')P!_ = [Fi(w) - F 2 (w) - F 3 (w) + F 4 (w)} P- P'_ . (96) 

In other words, F(v, v') actually may be taken to be a scalar, which we identify 
with the Isgur-Wise function, 

F(v,v') = £(w). (97) 

As an exercise, let us apply this formalism to the matrix elements for B — > 
(D, D*) iv. For a given matrix element, we pick out the part of the superfield 
M.(v) which is relevant. Hence we find 

(D(v')\ c^b \B(v)) = (M D (v')\ K^hl \M B (v) 

= -v^^iTr [ 7 5 p; 7 ^p + (- 7 5 )] £M 

= ^m D m B £(w) (v + v'f (98) 

(D*(v',e)\cn 5 b\B(v)) = (M* D (v',e)\h c v ,Y l l 5 h b v \M B (v)} 

= -^m D ,m B Tt yP^^P+i--/ 5 )] £(w) (99) 
= yJm D *m B £{w) \{w + 1)6*" -e*-(v + 
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(D*(v , ,e)\&fb\B{v)) = (M * D (v',e)\h c v ^h h v \M B (v)) 

= y/m D .m B £(w) ie^ ap tlv' a v p , (100) 



reproducing explicitly the relations fl23|) between the independent form factors 
hi{w). We can also derive the normalization condition at w = 1. Consider 
the matrix element of the b number current bj^b between B meson states. In 
QCD, the matrix element of this current is exactly normalized, 

(B(v)\ b^b \B{v)) = 2p B = 2m B v^ . (101) 

But in HQET, we have 

(B(v)\b^b\B{v)) = (M B (v)\h v ^h v \M B (v)) 
= m B ^(v ■ v)(v + v)^ 

= 2m B v IJ -£{l). (102) 
Hence the normalization condition at zero recoil, 

£(1) = 1, (103) 
follows directly from the conservation of the heavy quark number current. 

4-2 Power corrections to the matrix elements 

The matrix elements we have derived are computed in the strict limit mb. c — > 
oo. How are they affected by corrections of order 1/mb and l/m c ? There are 



two sources of 1/toq corrections in the effective theory: the corrections (64) 
to the heavy quark currents, and the corrections (|5^) to the Lagrangian. 

When radiative corrections are included, the expansion of the heavy quark 
current cTb in terms of HQET operators has a form which is somewhat more 
general than Eq. (|64[), 

cVb - a Q (a s ) h c v ,Th b v + ^2 h c v ,T«iD a h b v + 4^ K, (-iD ^V* h b v + .... 

/nib Zm c 

(104) 

The matrix elements of the power corrections are constrained by heavy quark 
symmetry in a manner completely analogous to the leading current. In terms 
of traces over the superfields, we have til 

(M'(v')\h v , T a iD a h v \M(v)) = -^M M M M - Tr ^m\v') T a M(v) G a (v, «')] , 

(105) 
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where G a (v, v') is another arbitrary 4x4 Dirac matrix. The matrix element 

(M'(v')\ h v , (-t < D a )T' 1 a h v \M(v)) (106) 

may also be written in terms of G a (v,v'), using charge conjugation. 

The 1 /rriQ corrections Ok and Oq to the Lagrangian contribute somewhat 
differently. In order to apply heavy quark symmetry, the matrix elements of 
the local currents, both leading and subleading, must be written in terms of 
the effective states \M(v)). However, these states are not eigenstates of the 
Hamiltonian, once Ok and Oq are included in the Lagrangian. Hence, for 
example, we must allow for the possibility that if an effective state \M(v)) is 
created at time t = -co, then Ok or Oq could act on the state before its decay 
at t = 0. This possibility is accounted for by including time-ordered products 
in which Ok or Oq is inserted along the incoming or outgoing heavy quark 
line. If we are keeping terms of order 1/mg, only one insertion qLOk or Oq 
needs to be included. The time-ordered products are of the formEil 

(D w (w)| cTb \B(v)) = 

•■■ + ^(M'iv'^ijdyT^ThlOK + O^lMiv)) 

+ ^(M'(v')\i J dyT{h c vl rhl,0 K + G }\M(v)), (107) 

where the ellipses include the current corrections computed earlier. The eval- 
uation of the matrix elements of the time-ordered products will lead to still 
more nonperturbative functions like F(v,v') and G a (v,v'). 

4-3 Corrections at zero recoil 

It is straightforward, but not very illuminating, to expand all of the new non- 
perturbative functions which arise at order 1/mg in terms of scalar form fac- 
tors. In the end, the corrections may be parameterized in terms of four func- 
tions of the velocity transfer w, and a single nonperturbative parameter A, 
all proportional to the mass scale Aqcdc-3 The new parameter has a simple 
interpretation as the "energy" of the light degrees of freedom, and is given by 

A = lim (m B - m b ) . (108) 

nit, — >oc 

Instead of a general treatment, however, we will consider the 1/mg cor- 
rections at the zero recoil point w = 1. This is clearly the most important 
case, because it is at this point that the nonperturbative matrix elements are 
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absolutely normalized in the heavy quark limit. What happens to this normal- 
ization condition when 1/toq corrections are included? 

Let us study the corrections to the current in detail. They are described by 
the nonperturbative function G a (v, v'). At v = v', G a (v,v) may be expanded 
as 

G a (v, v) = Gxv a + G 27q + G 3 v a tf + G 4l J . (109) 
But G a (v,v) is subject to the same constraint as F(v,v'), 

G a (v,v) = P^G a {v,v)P_ = (G 1 -G 2 -G 3 + G 4 )v a P_ = Gv a P_ , (110) 

and it, too, is equivalent to a Dirac scalar (the same is not true of the general 
function G a (v, v')). Now consider the matrix element where we take T a = v a . 
Then we have 



(M'(v)\h v ,iv-Dh v \M(v)) = - y/M M M M . Tr \m'(v) v a M(v) 



Gv a 



= - Gy/M M M M ' Tr M{v)M(y) . (Ill) 

But this matrix element vanishes by the classical equation of motion in the 
effective theory, 

v-Dh v (x) = 0. (112) 

Hence G = = G a (v, v). There are no lLajjQ. corrections from the current to 
the normalization condition at zero recoilr^r 1 ! 

The same is true of insertions of the corrections Ok and Oq to the La- 
grangian: their contribution vanishes at w = 1. To show this requires the 
imposition of the conservation of the b number current at order 1/m;,, much as 
we derived the normalization of the Isgur-Wise function at leading order. This 
part of the argument is_analogous to the classic nonrenormalization theorem 
of Ademollo and GattoE3 !_. 

In the end, we have the result known as Luke's Theoremt3 There are no 
corrections at zero recoil to the hadronic matrix elements responsible for the 
semileptonic decays B — > D Iv and B — > D* iD. The leading power corrections 
to the normalization of zero recoil matrix elements are only of order l/m 2 c . 
Given that Aqcd/^c ~ 30% and Aq CD /to^ ~ 10%, the implication is that 
the leading order predictions at w — 1 are considerably more accurate than one 
might have expected. In addition, away from zero recoil the 1 / m c corrections 
must be suppressed at least by (w — 1). 

On closer inspection, this result is more interesting for B — > D* lv than 
for B — > D lv. This is because the leading order matrix element for B — > D lv 
vanishes kinematically at zero recoil for a massless lepton in the final state. 
Hence, in this case the l/m c corrections are not suppressed as a fractional 
correction to the lowest order termed 
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4-4 Extraction of \V c b\ from B — > D* tv 

An immediate application of these results is the extraction of \V c b\ from the 
exclusive decay B — ► D* iv. This process is mediated by the weak operator 
Cfcc ((§)> whose matrix element factorizes as 



(D* tv\ O bc \B) = 



GpVcb 



(D* I gy (1 - j 5 )b \B) (£D\ £ ltl {l - r> |0) . (113) 



The leptonic matrix element may be computed perturbatively, while we treat 
the hadronic matrix element in the heavy quark expansion. The result is a 
differential decay rate of the forma 

7T" = ^\Vc b \ 2 (m B -m D ,) 2 m 3 D ,(w + l) 3 V^l 
aw 4»7r J 



1 



Aw m 2 B 



2wmi,mrj* + m 2 D , 



F 2 (w). 



(114) 



+ 1 {m B - m D ,) 2 

All of the HQET analysis goes into the factor F(w), which has an expansion 

F(w) = £(ui) + (radiative corrections) + (power corrections) . (115) 

We extract \V c b\ by studying the differential decay rate near w = 1, where 
the hadronic matrix elements are known. Of course, this requires extrapolation 
of the experimental data, since the rate vanishes kinematically at w = 1. For 
massless leptons, only the matrix element (D* \ cy^^b \B) of the axial current 
contributes at this point. The analysis of this quantity in the HQET yields an 
expansion of the form 



^(1) = VA 





+ + h/r. 

m c nib 



(116) 



The correction <5i/ m 2, which contains terms proportional to \/m 2 , 1/m 2 and 
l/m c mb, is intrinsically nonperturbatiue_.lt has been estimated from a variety 



(117) 



of models to be small and negativej---- 

S 1/m 2 w -0.055 ± 0.035 . 



Note that the model dependence in the result has been relegated to the esti- 
mation of the sub-subleading terms. The radiative correction tja has now been 
computed to two loopsp_'E_l 



t\a = 0.960 ±0.007. 



(118) 
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The result is a value for F(l) with errors at the level of 5%, 



F(l) = 0.91 ±0.04. (119) 

This is the theory error which the experimental determination of | T^t, | will 
inherit. It is dominated by the uncertainty in the nonperturbative corrections, 
and it is difficult to see how this can be improved much in the future. 

All that is left experimentally is to extrapolate the data to w = 1 and 
extract 

1 dr , s 

lim — . 120 

w-* 1 y/W — 1 dlV 



Once the kinematic factors in Eq. (114) have been included, this amounts to a 



direct measurement of the combination |V C &|F(1). Both CLEO and LEP have 
reported results for this quantityt^rcJ They have taken the slightly different 
value F(l) — 0.88 ± 0.05, and I have scaled up the theory error of the CLEO 
result to make it consistent with LEP. Then we have quite consistent results 
for Vet: 

CLEO: (39.4 ± 2.1 ± 2.0 ± 2.2) x 10 -3 
LEP average: (38.4 ± 1.1 ± 2.2 ± 2.2) x 10~ 3 . (121) 

This value of \V c b\ has almost no dependence on hadronic models. In contrast 
to model-based "measurements" , here the theoretical error is meaningful, in 
that it is based on a systematic expansion in small quantities. 



5 Inclusive B Decays 

An exclusive semileptonic B decay, such as B — > D £p, is one in which the final 
hadronic state is fully reconstructed. An inclusive decay, by contrast, is one in 
which only certain kinematic features, and perhaps the flavor, of the hadron 
are known. In this case, we need a theoretical analysis in which we sum over 
all possible hadronic final states allowed by the kinematics. Fortunately, this 
is possible within the structure of the HQET. 

As in the case of exclusive decays, the key theme is the separation of short 
distance physics, associated with the heavy quark, from long distance physics, 
associated with the light degrees of freedom. We will also rely on heavy quark 
spin and flavor symmetry. However, the new ingredient will be the idea of 
"parton-hadron duality" , which, as we will see, also relies on the heavy quark 
limit m b > A QC d- 
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5. 1 The inclusive decay B — > X c tv 



Let us consider the inclusive decay 



B(p B )^X e (p x )l(pi)v(p v ), 



(122) 



where all that is known about the state X c are its energy and momentum, 
and the fact that it contains a charm quark. This decay is mediated by the 
weak operator Ob c . It is easy to generalize our discussion to inclusive decays 
of other heavy quarks, such as b — > ulv and c — > s£v, by replacing 0\, c with 
the appropriate weak operator. 

The treatment of exclusive decays required both the b and c quarks to 
be heavy. For inclusive decays we can relax this condition on the c quark, 
requiring only mj, ^> Aqcd- What does the weak decay of the b, at time t = 0, 
look like to the light degrees of freedom? For t < 0, there is a heavy hadron 
composed of a point-like color source and light quarks and gluons. At t = 0, the 
point source disappears, releasing both its color and a large amount of energy 
into the hadronic environment. Eventually, for t > 0, this new collection of 
strongly interacting particles will materialize as a set of physical hadrons. The 
probability of this hadronization is unity; there is no interference between the 
hadronization process and the heavy quark decay. There are subleading effects 
in powers of AQCD/wf,, but they do not alter the probability of hadronization. 
Rather, they reflect the fact that the b quark is not exactly a static source of 
color: it has a small nonrelativistic kinetic energy and it carries a spin, both 
of which affect the kinematic properties of its decay. 

As in the case of exclusive decays, we will compute the inclusive semilep- 
tonic width . JXB. =r> A C ^P) as a double expansion in powers of a s (rrib) and 
A Q cD/m 6 E3BE3E3 The expansion in a s (nib) reflects the applicability of per- 
turbative QCD to the short distance part of the process. The heavy quark 
expansion will be continued to relative order 1/m^, as there is an analogue of 
Luke's Theorem which eliminates power corrections to the rate of order 1/mt,. 
These corrections will be written in terms of three nonperturbative parame- 
ters. The first, A, is defined in Eq. (|10S| ). It is essentially the mass of the 
light degrees of freedom in the heavy hadron, but we will see that it is plagued 
by an ambiguity of order Aqcd in the definition of the b quark mass. The 
other two parameters are the expectation values in the Sjpeson of the leading 




1 



(B\O k \B) 



2m B 
1 



(B\O g \B), 



(123) 



6ms 
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where we take the usual relativistic normalization of the states. Hence, Ai may 
be thought of roughly as the negative of the b quark kinetic energy, and A2 as 
the energy of its hyperfine interaction with the light degrees of freedom. 

Now let us outline the computation. The inclusive decay involves a sum 
over all possible final states, which is actually a sum over exclusive modes (such 
as D, D* , Dtt, . . .), followed by a phase space integral for each mode. We write 

T{B^X c lv)=Y J /d[P.S.] \ {X c lv\O bc \B)\ 2 . (124) 

There is an Optical Theorem for QCD, which follows from the analyticity of 
the scattering matrix as a function of the momenta of the asymptotic states. 
Its content is that a transition rate is proportional to the imaginary part of the 
forward scattering amplitude with two insertions of the transition operator, 

T(B^>X c £u) = -21m i j dx e tkx {B\T \o ] bc {x),O hc {^)] \B) = 2lmT. 

(125) 

In what follows, we will write the time-ordered product T{0\ c , Ob c } as a series 
of local operators, using the Operator Product Expansion. As we will see, the 
applicability of this expansion, and its computation in perturbation theory, 
will rest on the limit nib 3> Aqcd • We will then use this limit again to expand 
the matrix elements of these local operators in the HQET. 

The first step is to factorize the integration over the lepton momenta, 
which can be performed explicitly. Written as a product of currents, Obc takes 
the form 

Gj?Vcb 

where 



O bc = ^£j£ : J tli , (126) 



Jbc = *y"(i -7 5 )6 

J? = *y"(l - 7> . (127) 

Then T can be decomposed as an integral over the total momentum = 
Pe +Pu transferred to the leptons, 

T - \G 2 F \Vcb? J dqT^(q) V(g) . (128) 

Here the lepton tensor is 

L^M = J d[P-S.] (0| 4 \lv) (£9\ J tv |0> 

= ^~ (Qui* ~ 1 2 9^) , (129) 
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and the hadron tensor is 



T»»(q) =-ij dare** (B\ T {jg(x), J&(0)} \B) . (130) 

We will need the imaginary part, Im T^ v . Where is it nonvanishing? In 
quantum field theory, a scattering amplitude develops an imaginary part when 
there can be a real intermediate state, that is, the intermediate particles can 
all go on their mass shell. Whether this is possible, of course, depends on the 
kinematics of the external states. 

Ituthis case, there are two avenues for creating a physical intermediate 
statet3 The first is to act on the external state \B) with the transition current 
Jjf . The state which is created has no net b number and a single charm quark; 
the simplest possibility is the decay process b — ► c. The momentum of the 
intermediate state is px = Pb — Q] the condition that it could be on mass shell 
is simply 

Px = (PB - qf > m 2 D . (131) 
If we define scaled variables 

P B = niBV^ , q^ = q^/rriB, mo = mo/niB , (132) 
this condition becomes 

v -q< ^(l + q 2 ~m 2 D ) . (133) 

Another possibility is to act on \B) with the conjugate operator j£ c ]. This 
operation would produce an intermediate state with two b quarks and one c. 
For this state to be on shell, the momentum transfer has to satisfy 

Px = (PB + q) 2 > (2m B + m D ) 2 , (134) 

that is, 

vq> ~(3-f + 4m D +m 2 D ) . (135) 

The physical intermediate states are shown as cuts in the v ■ q plane in Fig. || 
Also shown is the contour corresponding to the phase space integration over 
the lepton momentum q. For physical (massless) leptons which are the product 
of a heavy quark decay, this integral runs over the top of the lower cut, for the 
range 

\/f + ie < v ■ q < - (l + q 2 - m 2 D ) + it . (136) 
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integration contour 

1 (extended) 

1 > K / 


v ■ q 




physical cut / 
unphysical cuts 



Figure 6: Analytic structure of T^ v in the complex v ■ q plane, for fixed, real q 1 . The 
integration contour is over the "physical cut" , corresponding to real decay into leptons. The 
unphysical cuts correspond to other processes. 

As indicated by the dotted line, we can continue this contour around the end 
of the cut and back along the bottom, to v-q = \ftp — ie. Since T^ u {v ■ q*) — 
—T^ u (v ■ q) for real q 2 , we compensate for extending the contour by dividing 
the new integral by two. 

We now encounter our central problem. The integral over v ■ q runs over 
physical intermediate hadron states, which are color neutral bound states of 
quarks and gluons. Hence the integrand depends intimately on the details of 
QCD at long distances, which is intrinsically nonperturbative. A perturbative 
calculation of T^, which is all we have at our disposal, would appear to be of 
no use. 

The solution is to deform the contour away from the cut, into the complex 
v ■ q plane, as shown in Fig. [7| Since the scale of momenta is setJay m&, 
the contour is now a distance of order Wf, away from the resonancesEHI Since 
nib Aqcd, it is reasonable to hope that a perturbative treatment in this 
region is valid. Essentially, we are saved because we do not need to know 
T^ v (q) for every value of q, just suitable integrals of T^. That we can use such 
arguments to compute perturbatively the average value of a hadronic quantity, 
where at each point the quantity depends on nonperturbative physics, is known 
as (global) parton-hadron duality. 

Parton-hadron duality has the status of being somewhat more than an 
assumption, since it is known to hold in QCD in the limit rrif, — > oo, but 
somewhat less than an approximation, since it is not known how to compute 
systematically the leading corrections to it. In any case, the limit mb ^> Aqcd 
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Figure 7: The deformation of the integration contour into the complex v ■ q plane. 

plays a crucial role here. By deforming the integration contour a distance 
of order m b away from the resonance regime, we find the correspondence in 
QCD of our earlier intuitive statement: the probability of the decay products 
materializing as physical hadrons is unity independent of the kinematics of the 
short distance process. The local redistribution of probability in phase space 
due to the presence of hadronic resonances is irrelevant to the total decay. 
Finally, we should note that since we do not have control over the corrections 
to local duality, it might work better in some processes than in others, for 
reasons that need not be apparent from within the calculation. Hence one 
must be particularly wary of drawing xlramatic conclusions from any surprising 
results of these inclusive calculationsE3 

Let us perform the operator product expansion at tree level, and for de- 
cay kinematics. The Feynman diagram is given in Fig. M, which yields the 
expression 

T v = fry"(l " 7 5 ) h-i + rnc _ 5 

(Pb - q) -mj.+ie 

We now write 

= m b v» + k" = 77i fc K + jfe") 
<f = q*/m b 
m c = m c /m b 

b(x) = e' mbV - x h v (x) + 0{l/m b ) , (138) 

and expand in powers of l/m b . Since the operator product expansion is in 
terms of the effective field h v , a factor of fc M corresponds to an insertion of the 
covariant derivative iD^. Operator ordering ambiguities are to be resolved by 
considering graphs with external gluon fields. 

As an example of the procedure, let us expand the propagator to order 
l/m b . (There are also corrections to the currents at this order, which are 
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Figure 8: The operator product expansion at tree level. 



included in a full calculation.) It is convenient to define the scaled hadronic 
invariant mass, 



s = {ni b v^ - q^ f/ml = 1 - 2v ■ q + q 2 
Then we find a contribution to T^ v of the form 



(139) 



T^ = — h v7 ^-$)y»(l- 7 5 ) 
mb 



1 



2k-q~k 2 



s — w? c + ie (s — m? c + ie) 2 



(140) 

From this expression we can read off the operators which appear in the operator 
product expansion. Since 

1 



Im 



raj + ie) 



= Tr(n-l)\(-l) n 5 ( - n - 1) (s-rh 2 c ), 



(141) 



we see that the effect of taking the imaginary part in each term is to put the 
charm quark on its mass shell. The leading term is a quark bilinear, 



— K'f(t-frf(l-'f)h v . 
nib 



(142) 



It is straightforward to compute its matrix element in the HQET using the 
trace formalism, 

(B|^7^^-^) 7 l/ (l-7 5 )/ lt ,|B) = 2m B (2i/V - g» v - v"q v - v v ef + g^v ■ q) . 

(143) 

The leading corrections to this expression are of order 1/m 2 ] there are no 
corrections of order l/rrib, by Luke's Theorem. Finally, we co ntract the tensor 
T^ v with L^ v and perform the phase space integration (128). In the end, the 
result is the same as we would have gotten directly by computing free quark 
decay, 

= G%\V cb \ 2 ml ( 

1927T 3 



(l - 8ml + 8m 6 c -ml- Ylm\ In m 2 c ) 



(144) 
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Of course, if we only intended to reproduce the free quark decay result, we 
would never have introduced so much new formalism. The value of the HQET 
framework is that it allows us to go beyond leading order and compute the 
next terms in the series in l/rri£. For example, consider the operators induced 
by the expansion of the propagator (14C). The correction of order l/m& comes 
from the operator 



h v j^-$)j u (l-j 5 )q-iDh v . 



m\ (s — m 2 + it) 2 
However, the matrix element of this operator is of the form 

(B\h v T a (v,q)iD a h v \B), 



(145) 



(146) 



which, as we have seen, vanishes by the classical equation of motion. (In 
writing Eq. (14C), we have already dropped terms explicitly proportional to 
v ■ k, for the same reason.) In fact, since all 1/mb corrections, from any source, 
have a single covariant derivative, they all vanish im the same way. This is 
the analogue of Luke's Theorem for inclusive decayso The correction of order 
1 1 m\ in Eq. ( pTp| ) is 



it)'- 



(147) 



The matrix element of this operator is related by the heavy quark symmetry 
to Ai, the expectation value of Ok- The full expansion of T^ v also induces 
operators with explicit factors of the gluon field, whose matrix elements are 
related to A2. 

We now present the result for the inclusive semileptonic decay rate, up 
to order 1/m 2 in the heavy quark expansion, and with the complete radiative 
correction of order a s . We also include that part of the two loop correction 
which is proportional to (3qCh 2 - Since /3q ~ 9, perhaps this term dominates the 
two loop result. In any case, it is interesting for other reasons, as we will see 
below. 

Let us first consider the decay B-^f^X u tv, for which the decay rate sim- 
plifies since m u = 0. We findEjl 



G%\Vub[ 

1927T 3 



1+ ^ 



25 



(2.98A, + C u ) 



2^2 
3 
2 



a s (m b ) 



(148) 



a s (m b )\ Ai - 9A 



+ 



2m 2 
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When wc include the charm mass, it is convenient to write the unknown quark 
masses in terms of the measured meson masses and the parameters of the 
HQET. In terms of the spin averaged mass Wlb = (tub + 3ms»)/4, we have 

m b = M B - A + + . . . , (149) 
2ni B 



and analogously for m c . We then nndBES£§E2 

.a s {m b ) 



T{B^X c tv)= ^m|x 0.369 
1927T d 



1 - 1.54- 



7T 



(1 .4 3/ 3 + C c ) \2*!*L) 2 - 1.65 A ( ! _ o.87^> 
7r / rriB V 77 



A 2 A A 

0.95^- - 3.18— i- + 0.02^- 

m B m B m B 



(150) 



All the coefficients which appear in this expression are known functions of 
mo/ruB, and are evaluated at the physical point mo/ms = 0.372. In both 
B — > X u tv and B — > X c £p, the power corrections proportional to Ai and A2 
are numerically small, at the level of a few percent. 

5.2 Renormalons and the pole mass 

The inclusive decay rate depends on the heavy quark mass mb, either explicitly, 



as in Eq. ( |148| ), or implicitly through A, as in Eq. (150). At tree level, mb is 
just the coefficient of the bb term in the QCD Lagrangian, but beyond that we 
are faced with the question of what exactly we mean by mb- Should we take an 
MS mass, such as mjfmj)? Or should we take the pole mass TO^ ole , or maybe 
some other quantity? The various prescriptions for mb can vary by hundreds of 
MeV, and, since the total rate is proportional to m\ , the question is of practical 
importance if we hope to make accurate phenomenological predictions. 

At a fixed order in QCD perturbation theory, the answer is clear. The 
heavy quark masses which appear come from poles in quark propagators, so 
we should take m^ olc (and mj? ole ). This is also the prescription for the mass 
which cancels out the on-shell part of the heavy quark field in the construction 
of £hqet- Hence the difference of heavy quark pole masses is known quite 
well, 

< ol ° - < olc = (m B - A + -^i- + ..) - (m D -A 1 A| 



2ms J \ 2mu 
3.34GeV + 0(A 2 QCD /m 2 3 ). (151) 
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Figure 9: The radiative corrections to m.^° lc of order a s (mj,)™" 1 " 1 /?^ . 



Since T(B — > X c £i>) depends approximately as m%(mb — m c ) 3 , the uncertainty 
due to quark mass dependence is reduced. 

The problem, of course, is that there is no sensible nonperturbative defi- 
nition of m^ olc , since due to confinement there is no actual pole in the quark 
propagator. Hence a direct experimental determination of a value for mj°'° to 



insert into the theoretical expressions ( |148| ) and (150) is not possible. How, 
then, can we do phenomenology? 

One approach would be to define rn^ olc to be the pole mass as computed 
in perturbation theory, truncate at some order, and then estimate the theo- 
retical error from the uncomputed higher order terms. However, it turns out 
that even within perturbation theory the concept of a quark pole mass is am- 
biguous. Consider a particular class of diagrams which contribute to m% olc , 
shown in Fig. [|. The perturbation theory is developed as an expansion in the 
small parameter a s (mb), so we hope that it will be well behaved. Each of the 
bubbles represents an insertion of the gluon self-energy, which is proportional 
at lowest order to a s (mb)^o- Of course, the infinite sum of the graphs in Fig.|| 
can be absorbed into the one loop graph, with a compensating change in the 
coupling from a s (rrib) to a s (q), where q is the loop momentum. The result is 
an expansion for m? olc of the form 



pole / \ 

ml =mb(mb) 



1 + aia s + (a 2 Po + b 2 )a 2 s + (a 3 /3g + b 3 /3 + c 3 )a 3 s 



, (152) 



where a s = a s (m,b). The graphs in Fig. |^ contribute the terms proportional 
to a™ +1 /3g . Since 0o « 9 these terms are "intrinsically" larger than ones with 
fewer powers of /?o, and we might hope that their sum approximates the full 
series. However, it is important to realize that the only limit of QCD in which 
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such terms actually dominate is that of large number of light quark flavors, 
in which case the sign of /3q is opposite to that of QCD. Although this is a 
physical limit of an abelian theory, we are certainly not close to that limit here. 
The ansatz of keeping only the terms proportional to (nib) (3q is known as 
"naive nonabelianization" (NNA)c3 

What is most interesting about the series of terms shown in Fig. |^, which 
takes the form ^a n a™/3o _1 , is that it does not converge. Already in the 
graphs kept in the NNA ansatz, we are sensitive to the fact that QCD is an 
asymptotic, rather than a convergent expansion. For large n the coefficients 
a n diverge as n\, much stronger than any convergence due to the powers a™. 
The series can only be made meaningful if this divergence is subtracted. As 
with many subtraction prescriptions, there is a residual finite ambiguity.^] This 
ambiguity, known .as. an "infrared renormalon" , leads to an ambiguity in the 
pole mass of orderE3'E3c3 

6ml olc ~ 100 McV. (153) 



By the definition (108), A also inherits this ambiguity. 



The expressions (148) and ( |150] ) are plagued by two problems. The first is 
the renormalon ambiguity in c and A. The second is that the perturbative 
expansion for the rate T is itself divergent, and also has an infrared renormalon. 
In the expansion 

T = T ^^a^a™(TO b )/3Q _1 + (power corrections) , (154) 

the coefficients a' n also diverge as n\. However, it turns out that these two 
problems actually cure each other, because, the infrared renormalons in m^ ole 
and in the perturbation series for V cancelu-3'c!3 We can exploit this cancelation 
to improve the predictive power of the theoretical computation of the rate. 
Without this improvement, the infrared renormalons render the expressions 



(148) and (15C) of dubious phenomenological utility. 



The most reliable approach, theoretically, is to eliminate rn^ olc or A ex- 
plicitly from the rate by computing and measuring another quantity which 
also depends on it. For example, let us consider the charmless decay rate 
T(B — > X u tv) and the average invariant mass (sh) of the hadrons produced 
in the decay. Each of these expressions suffers from a poorly behaved per- 
turbation series in the NNA approximation. Ignoring terms of relative order 
1/m^ and writing the rate in terms of A instead of m^ olc , we find to five loop 



d In a formal treatment, this ambiguity arises from a choice of contour in the Borel plane. 
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orderC 



_ G 2 F \V ub \ 2 5 
1 ~ 192^ m B 



Gp\V u b\ 2 5 

£ n 

192tt 3 s 



1 - 2.41— - 2.98 

7T 



(f) A-44,(i)'/« 

»(7)''S-'»(7)'fr-4 



1 - 0.061 - 0.120 - 0.107 - 0.111 - 0.136 - 
- 5A/m B + . . 



(155) 



for a s (mb) = 0.21 and fto = 9. As we see, not only does the perturbation series 
fail to converge, it does not even have an apparent smallest term, where one 
should truncate to minimize the, error of the asymptotic series. The series for 
(sh) exhibits a similar behavioral 



).20^ + 0.35 (^) 2 A) + 0.64 fa + 1.29 (^)* /3 3 



2.95 



(-) 5 A « + ... + -LA + 

V 7r / 10 ms 

= m| [0.0135 + 0.0141 + 0.0156 + 0.0189 + 0.0261 
-7A/10m B + ...] . 



(156) 



However, the situation improves dramatically if we eliminate A and write T 
directly in terms of (sjj), 



G 2 F \V u b\ 2 _ 5 
192tt 3 B 



1 - 7.14 



- 0.064 - 0.020 - 0.0002 - 0.022 - 0.047 - 



(157) 



By truncating this series at its smallest term, 0.0002, we obtain a new expres- 
sion in which the theoretical errors are under control. The price is that we 
must now measure a second quantity, (sh), in the same decay. 

In principle, the same procedure works for decays to charmcJ In practice, 
it is best to combine a number of determinations of A. This has been done by 
CLEOpH which has performed a comparison of measurements of the moments 
(Eg), (E 2 ), (s H — rn 2 D ) and ((sh —m 2 D ) 2 ). The results, reproduced in Fig. ITol, 
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A (GeV) 

Figure 10: The CLEO determination of A and Ai from B — » X c £u. 

are somewhat disappointing, in that one does not obtain a very consistent 
determination of A and Ai. Perhaps this is just a fluctuation, or perhaps this 
is a sign that parton-hadron duality is failing in these days. 

An alternative approach is to express the width Xjiu terms of the running 
mass m&(m&) instead of another inclusive observableOcH Since the MS mass is 
a short distance quantity, this also eliminates the infrared renormalon, which 
is associated with long distance physics. However, from a phenomenological 
point of view, it raises the question of how the running mass is to be deter- 
mined from experiment. Possibilities include quarkonium spectroscopy, QCD 
sum rules, and lattice calculations, but in all of these cases it is important to 
determine reliably the accuracy of the method, and how to deal with renor- 
malon ambiguities in a manner that is consistent with their treatment in the 
calculation of T. Nevertheless, such an approach, particularly one based on 
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analyzing quarkonium and production near threshold, should eventually prove 
fruitfully 



5.3 Phenomenology of V c b and V u b 

Despite this ambiguous situation, groups have presented extractions of V c b 
based on inclusive semileptonic b decays by simply inserting "reasonable" val- 
ues for A and Aj_. Such an approach clearly has its dangers! At any rate, the 
quoted result isc3 



The lack of controversy about this procedure is no doubt due in part to the fact 
that this number is quite consistent with that determined from the analysis of 
exclusive decays. 

There are additional, much more interesting, problems with extracting 
\V u b\ from the inclusive decay B — > X U £D. They arise from the problem that 
the process B — > X c £u presents an overwhelming background to B — > X u lv. 
The only way to avoid this background is to restrict oneself to a corner of phase 
space in which charmed final states are kinematically inaccessible. Existing ex- 
perimental analyses isolate the charmless decays by imposing the requirement 
Ef > {m 2 B — mj D )/2mB or s H < m 2 D . Unfortunately, the OEE_can be shown 
to break down in these restricted corners of the phase spaceCJtiEl 

I do not have space here to explore this issue in much detail, but it is easy 
to appreciate the essence of the problem. For massless final states, the OPE 
is an expansion in powers of the light quark propagator, l/m&(l — v ■ q + q 2 ). 
Over most of the final state phase space, the denominator is of order nib and 
the OPE is well behaved. But there exist configurations for which both the 
denominator vanishes and the operator matrix elements which appear in the 
OPE are nonzero. It turns out that the dangerous region is when v ■ q — > \ 
and q 2 — > 0. The precise form of the divergence depends on the kinematic 
distribution being studied. The general form, however, is universal. In this 
"endpoint region" , let y be a scaled variable such that y — > 1 at the kinematic 
endpoint. For example, for dT/dEe, we take y = 2E(/mb, and for dT/dsH wc 
take y = sg/Arai,. Then near y = 1, the OPE takes the general form 



(B(v)\ h v iD m . . . zZV hy \B(v))/2m B = A n . . .v^ + . . . . (160) 



\V cb \ = (40.0 ± 0.4 ± 2.4) x 10~ 3 . 



(158) 




(159) 



where 



c n are coefficients of order one, and the A n are moments defined by 
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The ellipses represent terms involving factors of the metric tensor , which 
are subleading. Since the A n are associated with totally symmetric combi- 
nations of the covariant derivatives, they may be interpreted roughly as the 
moments of the heavy quark momentum. Note that as defined, Aq = 1, A\ = 
and A2 — Ai. 

The divergences in Eq. ( |l59| ) can be controlled only if one integrates over a 
large enough region near the endpoint, 1 — S < y < 1. For S ~ A/ mj,, one finds 
a series for F which does not converge, since the individual terms are of order 
A n /A n ~ 1. This situation reflects a dependence of the shape of dT/dy on the 
entire b quark momentum distribution in the B meson. Since, for example, the 
window 2.3 GeV < Eg < 2.6 GeV corresponds to 6 ~ 0.1, this problem pollutes 
the extraction of \V u b\ from dT/dEg. One is forced to introduce a model for 
the b wavefunction, with the attendant uncontrolled theoretical uncertainties. 
The same is true, it turns out, for dT/dsH- The current "best" measurement 
of V u b from LEP, based on an inclusive analysis, isEa 

\V ub \ = [4.05i°-^(stat.)i°-gl°-g(sys.) ± 0.02fo) ± 0.16(HQE)] x 10~ 3 , 

(161) 

or approximately IK^/Vy = 0.104+qq]J. While these analyses are experimen- 
tally very sophisticated, they rely intensively on a two-parameter model of the 
6 quark wavefunction. Essentially, in such a parameterization all moments of 
the b momentum distribution are correlated with the first two nonzero ones, a 
constraint which is unphysical. Even if the two parameters are varied within 
"reasonable" ranges, it is doubtful that such a restrictive choice of model cap- 
tures reliably the true uncertainty in \V u b\ from our ignorance of the structure 
of the B meson. While the central value which is obtained in these analyses is 
reasonable, the realistic theoretical error which should be assigned is not yet 
well understood. ._. 
A recent analysis by CLEO of the exclusive decay B — > plv yields 

IK&I = [3.25 ± 0.14(stat.)l°;^( S yst.) ± 0.55(theory)] x 10" 3 , (162) 

or approximately \V u b/V c b\ = 0.083tooi6i essentially consistent with the LEP 
result. In this case the reliance on models is quite explicit, since one needs the 
hadronic form factor (p\ -07^(1 — r y 5 )b \B) over the range of momentum transfer 
to the leptons. The CLEO measurement relies on models based on QCD sum 
rules, which have uncertainties which are hard to quantify. Hence, just as 
in the case of the LEP measurement, the quoted errors should not be taken 
terribly seriously. All of the current constraints are consistent with |Kjf,/V^| = 
0.090 ± 0.025, where I strongly prefer this more conservative estimate of the 
theoretical errors. The problem lies not in the experimental analyses, but in 
our insufficient understanding of hadron dynamics. 
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Recently it has been pointed outE.3 that the one may alternatively reject 
the charm background by studying the distribution dT/dq 2 and restricting 
oneself to q 2 > (rn,B — itt-d) 2 ■ Not only does this cut eliminate charmed final 
states, but it also avoids the troublesome region near q 2 — 0. Hence such 
a determination would not be polluted by the divergence described above. 
Whether the neutrino reconstruction algorithms of the B Factories will be up 
to the task of this measurement yet remains to be seen. 



6 Concluding Remarks 

Unfortunately, we have had time in these lectures only to introduce a very 
few of the many applications of heavy quark symmetry and the HQET to the 
physics of heavy hadrons. Since its development less than ten years ago, it has 
become one of the basic tools of QCD phenomenology. Much of the popular- 
ity and utility of the HQET certainly come from its essential simplicity. The 
elementary observation that the physics of heavy hadrons can be divided into 
interactions characterized by short and long distances gives us immediately a 
clear and compelling intuition for the properties of heavy-light systems. The 
straightforward manipulations which lead to the HQET then allow this intu- 
ition to form the basis for a new systematic expansion of QCD. The deeper 
understanding of heavy hadrons which we thereby obtain is increasingly im- 
portant as the B Factory Era begins. 
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